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Abstract 

The forces and moments at the rotor hub of an Aerospatiale SA349/2 helicopter were investigated. 
The study included three main topics. First, measured hub forces and moments for a range of level flight 
conditions (n;=0.14 to 0.37) were compared with predictions from a comprehensive rotorcraft analysis to 
examine the influence of the wake model on the correlations. Second, the effect of changing the blade 
mass distribution and blade chordwise center of gravity location on the 3/rev nonrotating frame hub 
loads was studied for a high-speed flight condition (H;=0.37). Third, the use of higher harmonic control 
to reduce nonrotating frame 3/rev hub shear forces was investigated. The last two topics were 
theoretical studies only. 
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chordwise blade aerodynamic center location 
number of blades 
blade chord, m 

rotor blade coordinate system; c-parallel to blade chord, in hub-plane, parallel to 
blade radius, in hub-plane, n-normal to hub-plane 
controller gain (response feedback) 
chordwise blade center of gravity location 
nonrotating frame x coefficient of hub force 
nonrotating frame y coefficient of hub force 
nonrotating frame z coefficient of hub force 

rotor thrust coefficient divided by rotor solidity, rotor thrust/TcR^p(QR)^a 
local blade drag, N 

single-blade hub force in the c-direction, N 

aerodynamic component of fc , N 

inertial component of fc , N 

D/q of the rotor hub, m^ 
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aerodynamic component of fn, N 
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component of fn due to first elastic torsion (k=l) modal component, N 
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total nonrotating frame hub force in the Z-direction (Fz=ftz)^ N 

acceleration of gravity, m/s^ 

structural damping coefficient 

characteristic blade inertia, m^kg 

blade modal inertia 

quadratic performance function for higher harmonic control, eq. 12 

distance from hub plane to shaft strain gages at 'a', m 

distance from hub plane to shaft strain gages at V, m 

local blade lift, N 

sectional blade mass, kg/m 

measured shaft bending moment at ’a’, Nm 

measured shaft bending moment at 1?’, Nm 

rotor hub mass (total rotor mass - total blade mass), kg 

aerodynamic moment (forcing function of blade flap/lag equation of motion), Nm 
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total rotating frame hub moment in the x-direction^ Nm 

total rotating frame hub moment in the y-direction, Nm 

total rotating frame hub moment in the z-direction, Nm 

total nonrotating frame hub moment in the X-direction, Nm 

total nonrotating frame hub moment in the Y-direction, Nm 

total nonrotating frame hub moment in the Z-direction Nm 

generalized blade torsion degree of freedom 

generalized blade bending degree of freedom 

rotor radius = 5.25 m 

blade radial location, m 

blade lag inertial constant 

blade flap inertial constant 

transfer-function matrix 

component of local blade velocity along radius (normal to the c-n plane), m/s 

local blade velocity, m/s 

flight velocity, m/s 

rotating frame shaft axis system 

nonrotating frame shaft axis system 

vector of harmonics of response variables, eq. 10 


a blade angle of attack, deg 

Os shaft angle of attack, deg 

Pic longitudinal flap angle, deg 

T|i lag bending mode shape 

T|f flap bending mode shape 

(]) induced angle of the resultant air velocity in the c-n plane, deg 

y lock number, pacR^/Ib 

\i advance ratio = Vq/QR 

V blade bending modal frequency, 1 /rev 

'P blade azimuth angle, measured from tail boom toward advancing side of rotor disk, deg 

p air density, kg/ m3 

a solidity, bc/icR 

00 collective pitch, deg 

0c lateral cyclic pitch, deg 

vector of harmonics of control variables (fixed-frame), deg 
0e elastic pitch, deg 

0HfIC HHC control vector (rotating frame), deg 

0p pilot control input, deg 

0s longitudinal cyclic pitch, deg 

^ torsion mode shape 

Q rotor rotation rate, rad/sec 


Subscripts 

c cosine component 

i harmonic number 

k mode number 

n harmonic number 

s sine component 

1. Introduction 

Reducing helicopter vibration is desirable for many reasons, including increased crew and 
passenger comfort, and reduced wear on the aircraft An understanding of the sources which generate 
vibration can aid in the development of rotorcraft with lower oscillatory loads. Examining helicopter 
forces and moments at the rotor hub, both experimentally and theoretically, is an important step 
toward understanding vibration. With this in mind, three aspects of hub loads were studied. The first 
topic of the paper compares predicted and measured hub loads in both the rotating and nonrotating 
frame. Predictions were made using both a uniform inflow and a nonuniform inflow (prescribed wake) 
model to determine the effect of the wake model on the correlations. The second topic investigates how 
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variations in the blade mass distribution and chordwise eg location influenced the SA349/2 oscillatory 
hub loads. The topic of altering blade design to reduce hub loads has been investigated both 
experimentally and analytically. References 1-8 address the effects of torsional stiffness, tip sweep, 
camber, chordwise eg offset (relative to the ac axis), blade mass distribution, twist, blade bending 
stiffness, and tab deflections on hub vibration and blade vibratory loads of four-bladed articulat^ 
rotors. This paper concentrates on predicting and studying the oscillatory hub loads of a thiee-bladed 
helicopter. The third topic investigates using higher harmonic control (HHC) of the rotor to minimize 
the hub shear forces. The last two topics were strictly theoretical investigations. 

This effort represents the second phase of a cooperative program between the National 
Aeronautics and Space Administration (NASA) and the French Ministry of Defense. The purpose of 
this program is to learn more about the aerodynamic and dynamic behavior of helicopters using 
SA349/2 flight test data and current helicopter analyses developed in the United States and in France. 
The first phase began with the 1984 flight test of the SA349/2 helicopter. Rotor blade structural and 
aerodynamic loads were measured for various thrust levels and advance ratios. Reference 9 presents 
the flight test data, and includes a description of the data acquisition and reduction procedures. 
References 10 and 11 present correlations between analyses and the aerodynamic and structural data. 
The effects of wake geometry, dynamic stall model, transonic aerodynamics, and different methods of 
calculating the blade elastic response were investigated. Reference 12 concentrated on the effect of 
individual blade modes on blade aerodynamic and structural loads. Reference 13 validated a new 
method of computing the forced response of an elastic blade using the measured aerodynamic loads 

The SA349/2 aircraft was flown again in 1987, thus initiating the second phase of this joint 
program. For this flight test, the rotor mast was instrumented to measure hub forces and moments. More 
detailed blade bending moment data were also obtained. Several 1984 flight conditions were 
duplicated during the 1987 flight tests by matching ^and Cj/a For these duplicated conditions data 
for blade aerodynamic and structural loads, hub loads, and rotor flapping angle and shaft angle are 
available. Collectively, results from the two flight tests provide a comprehensive data base. Data 
from the 1987 flight test are documented in ref. 14. 

2. SA349/2 Flight Test Data 

The SA349/2 aircraft, fully described in ref. 9, has an advanced geometry main rotor and three 
Grande Vitesse blades with rectangular tips. During the 1987 flight test, the main rotor was 
instrumented with 39 hub and blade strain gages, and six shaft strain gages. Descriptions of the 
instrumentation, the data acquisition and the data reduction procedures are published along with the 
data in ref. 14. Table 1 lists the basic SA349/2 aircraft characteristics. The six force and moment 
components at the hub were obtained from strain-gage measurements on the rotor shaft. Four of the six 
shaft strain gages measured the shaft bending moments, m^ and m\^, at distances la and lb below the 
rotor (fig. 1). The measurements of ma and mb were sufBcient to derive the total in-plane rotating frame 
forces and moments at the hub (ftx, fty, nttx, nity). This was possible since the x- and y-components of the 
forces and moments, although 90 deg. out of phase, experienced the same loading conditions around the 
eizimuth and thus were of equal magnitude. The in-plane rotating frame forces and moments at the hub 
were calculated assuming the shaft was a beam isolated from other external forces and moments. In 
deriving the rotating frame hub loads, then, the x-component was computed from the measured shaft 
loads, and the y-component was obtained by introducing a 90 deg. phase shift. The nonrotating frame 
forces and moments were derived as outlined in sections 3.1 and 3.2. The two remaining shaft gages 
measured rotor thrust and rotor torque, both of which are the same in the rotating and fixed frame. 

3. Analytical Approach 

Three topics related to the SA349/2 hub loads were studied using a comprehensive rotorcraft 
analysis. The first was a correlation of SA349/2 hub loads flight test data with predicted hub loads, to 
assess the effect of the wake model on the predictions. Second, variations in the hub forces resulting 
from alterations in the blade radial mass distribution and chordwise eg location were examined. 

Finally, the effect of using HHC to reduce predicted SA349/2 hub forces was investigated.The last two 
topics were theoretical investigations. 

The SA349/2 hub loads were studied using a comprehensive analytical model of rotorcraft 
aerodynamics and dynamics (CAMRAD). CAMRAD (refs. 15-16) is an analysis used to examine a 
variety of two-rotor helicopter designs. For this study, predictions of the rotor hub forces were obtained 
by trimming the CAMRAD model of the SA349/2 to the aircraft forces and moments. Five blade 
bending modes; rigid flap and lag , and three coupled flap/lag elastic bending modes, and two elastic 
torsion modes were included in the rotor analysis. Control system stiffness, essentially a rigid torsion 
mode, was included as well. 
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3 1 De rivation of Total Rotating and Nonrotnting Frame Hub Forces and Moments. 

The rotating and nonrotating frame hub forces and moments are shown in fig. 1. Tlie components ftx, 
fty, nitx, and mty represent the combined in-plane forces and moments for all three blades in the rotating 
frame, shaft axis system (x, y, z). As mentioned in section 2, the in-plane rotating frame forces, f^x and 
fty , have equal amplitudes. The same is true for the in-plane moments m^x and mty. The components 
Fx/ Fy/ Mx/ and My are the corresponding in-plane loads in the nonrotating frame, shaft-axis system 
(X, Y, Z). Fz(=ftz) and Mz(=mtz) are the out-of-plane loads in the nonrotating frame, shaft-axis 
system. 

The total fixed-system hub forces were obtained from the total rotating hub forces, which in turn 
were derived from the blade root forces, as shown in fig. 2. Figure 2a shows the single-blade hub force 
components, which are obtained by integrating the local forces along the blade. The three force 
components acting at the root of each blade are fc the hub-plane force in the chordwise direction, fr, the 
hub-plane force in the radial direction, and fn, the blade force normal to the hub-plane. Since these 
force components are functions of time, they can be represented for the blade as a sum of i harmonics; 

fc('f^) = Z fed cos(iH^) + fcsisin(i4*,) 


fr(^) = I frd cos(i'f^) + f„i sinfiTj) 


( 1 ) 


fn(^P = S fnd cos(iT^) + fnsi sin(i'Pj), 

where the subscripts ’ci’ and si’ refer to the amplitudes of the cosine and sine components of the i^^ 
harmonic. Rgure 2b shows the force components for all three blades, where blade one is at deg 

(aligned with the tail boom), blade two is at 4^=120 deg and blade three is at 'F3=240 deg. The total 
rotating frame hub force components in the x, y, z axis system are then obtained by summing the single- 
blade components as follows; 

ftxW = W'Fi) - frf'Fi+nOsinSO - fr('Fi+240)sin30 

- fc(H'i+120)cos30 + fc('Fi+240)cos30 

fty(40 = -fc('Fi) + fcf'Fi+nWsinSO + fc(Ti+240)sin30 (2) 

- fr('Fi+120)cos30 + ^4^1+240)00830 

ftz(40= fn('Fi) + fn('Pi+120) + fn(4'i+240). 

Substituting eqs. 1 into eqs. 2 and transforming from the single-blade coordinate system (c, r, n) to the 
total rotating frame coordinate system (x, y, z) yields equations for f^x/ fty/ and ftz in harmonic form (ref. 
14). Interestingly, because of cancellations due to rotor symmetry and relative blade phase, only the 
nb+1 harmonics of fix and fjy are non-zero, and only the nb harmonics of ftz exist in the rotating frame. 
When transforming to the fixed system, the nb-1 and nb+1 harmonics of the in-plane forces in the 
rotating frame combine to produce only nb harmonics of Fx and Fy. These fixed system forces are 
obtained through the following coordinate transformation, pictured in fig. 2c; 

Fx(40= ftxcos(40 + ftysin(49 

Fy(40 = -ftx sin(40 + fty cos(40 (3) 


Fz(40= ftzW. 

The results section presents correlations of the in-plane loads with forces and moments in both the 
rotating and nonrotating frames. Because the nb-1 and nb+1 harmonics in the rotating frame combine to 
produce the nb harmonics in the fixed-frame, the relationship between the rotating and nonrotating 
frame forces is not straightforward. In addition, aircraft vibration is generated by the fixed frame 
forces but these forces originate in the rotating frame. Thus it is instructive to examine results in both 
the rotating and nonrotating frames. 
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3.2 Derivation of Hub Force Aerodynamic and Inertial Components. 

To determine the mechanisms through which the hub forces were altered by changes in blade mass 
distribution and blade eg location, the components of the single-blade forces, fc, fr, and fn, which make 
up the total forces as shown in fig. 2, were examined. The single-blade in-plane forces can be 
decomposed into aerodynamic and inertial components using the following relationships (ref. 17, pp. 64, 
100 - 102 ), 


fc — fca + td 

= { r (Lsini|) + Dcos<t»)dr/R} + {[l/(2aR/oll2l{,)][S SqU^(qi(-qk)]} (4) 

and •' 

fra + frf 

= {J[Dur/U-Sqk(nfkfna-illkfca)Wr/R) + {[l/(2aR/a£)?Ib)](S 2Sqikqk)) (5) 

where fca and f^ are the aerodynamic components of fc and ff, and fei and fn are the inertial 
components. The aerodynamic chordwise force, fca/ is a function of the c-components of blade lift (L) and 
drag (D); fra is a function of blade drag, the r-component of the local blade velocity (ur/U), the modal 
displacements (qk^lfk and qk^llk/ where qk are the nondimensional blade degrees of freedom,‘nfk and qik 
are the flap and lag components of the bending mode shapes), and the c- and n-components of the single- 
blade aerodynamic forces (fca and fna)- Tl'® components fd and frf are, neglecting hub and shaft motion, 
simply functions of blade lag inertial constants, Sq|k/ and nondimensionalized blade bending degrees of 
freedom, qk- The subscript k refers to the bending modal component. The individual inertial and 
aerodynamic contributions to f^x obtained by combining eqs. 4 and 5 according to the 

relationship given in eq. 2. 

The 3/rev out-of-plane forces can be similarly broken down into aerodynamic and inertial 
components using the equation (ref. 17, pp. 65, 100), 

fn= fna + frf 

= { J (Lcos(|) - Dsin(t.)dr/R) + {[l/(2aR/ofl2lb)](I Sqfkqk)) , (6) 

where f^a/ the aerodynamic component of fn, is a function of the components of lift and drag; and fni, the 
inertial component of fn/ is a function of the blade flap inertial constants, Sqfk/ ^nd the 
nondimensionalized blade degrees of freedom, qk* The total out-of-plane forces for the inertial and 
aerodynamic components are then derived according to the relationship in eq. 2 for f^z. 


The inertial components of the forces, fei, fn, and fni, are comprised of the inertial constants and 
the nondimensionalized blade degrees of freedom. The inertial constants are given by (ref. 17, p. 69) 

Sqlk = (l/Ib)^‘nikoidr and Sqflc = (l/Ib) nfkn^dr (7) 

where Ib is the characteristic blade moment of inertia, T|ik and *nfk are the lag and flap components of the 
bending mode shapes, and m is the sectional blade mass. The nondimensional blade degrees of freedom 
are derived from the following coupled flap/lag equation of motion (ref. 17, p 62), 


lqk(qk + gs^kqk ^k^k) + (nonlinear terms, torsion coupling,...) 

= 7(l/ac)Mqk(aero) + constant, (8) 

where Iqk = (1/Ib) t|k^nidr/R is the blade modal inertia. Thus the nondimensionalized blade degrees 
of freedom are functions of bending modal frequencies, Vk/ modal inertias, Iqk/ and aerod)mamic forcing 
functions, (l/ac)Mqk(aero)- The aerodynamic forcing functions can be described by the equation (ref. 17, 
p. 67), 

(l/ac)Mqk(acro) = J (Vac)(T|fkfna-nikfca)dr. (9) 

The inertial components of eqs. 4-6 suggest that any changes in the inertial hub forces caused by 
altering the blade properties will arise through changes in the bending modal characteristics (vk, Iqk/ 
Tifk, Tiik) and/or the aerodynamic forcing function ((l/ac)Mqk(aeix>))* The blade bending modal response 
depends partly on the blade mass distribution, so altering the blade mass distribution is expected to 
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affect the hub loads through changes in the bending modal response. In contrast, eq. 9 indicates that 
the aerodynamic forcing function depends not only on the bending modal characteristics, but also on the 
blade lift through fna and fca» Likewise, the aerodynamic components of eqs. 4-6 are all influenced by 
the blade lift. The lift is a function of the blade angle of attack, a, which depends on the blade pitch 
(both rigid and elastic) and the local blade velocity (ref. 12). Changes in the aerodynamic components 
of the hub forces caused by altering the blade properties can be traced to changes in the elastic 
component of blade pitch, 0e. The elastic blade pitch is a function of the elastic torsion modal response 
of the blade. Since the chordwise eg location affects this dynamic torsion response, shifting the eg 
location is expected to affect the hub loads through changes in the blade angle of attack. Section 4.2 
explores these speculations. 

3.3 Higher Harmonic Control Algorithm 

The effectiveness of using HHC to reduce the SA349/2 hub loads was investigated with 
CAMRAD/JA, a modified version of G\MRAD developed at Johnson Aeronautics. The objective was to 
use 3/rev swashplale control to minimize the 3 /rev hub shears. The measurement vector was given by, 

z = I(Cfx/o)3c (^FX/^bs (Cfy/<^)3c (Cfy/<^)3s (Cpz/<^bc (Cfz/<^)3s1/ HO) 

where CpX/ Cpy, and Cpz the nonrotating frame hub force coefficients relative to the shaft axes, a is 
the rotor solidity, and the subscripts *3c' and *3s* refer to the 3 /rev cosine and sine components. The 
control vector was given by, 

®con ^ [(^)3c (®o)3s (®c)3c (®c)3s (®s)3c (®s)3sl/ (11) 

where Gq, 0c, and 0s are the nonrotating frame collective, lateral cyclic and longitudinal cyclic pitch 
controls. 

A controller was introduced between the trim and circulation /motion iteration sections of the 
solution procedure. The controller was similar to the self-tuning regulators used for HHC of helicopters 
(refs. 18 and 19). However, instead of applying the HHC algorithm to a dynamic system (as would be 
the case in flight), a Newton-Raphson algorithm was applied to the aircraft at a spedtied operating 
condition. As a result, the controller became just an algorithm to Bnd the optimum solution to a 
nonlinear algebraic problem. 

The controller uses the quadratic performance function, 

J = z^z + O.OOO2A0con ^A0oon / 02) 

where A0con , the change in controls from one iteration to the next, was introduced to improve 
convergence. A locally linearized model of the helicopter, Az=TA0con / ''^^s assumed (ref. 18). T was 
the transfer-function matrix. Minimizing the dimensionless performance function, J, led to the linear 
feedback equation A0con =-Cz, where C was the controller gain. Convergence of the controller was 

defined in terms of the difference between successive iterations of J, where AJ < l.OE-8 indicated 
convergence. At the beginning of each Newton-Raphson trim iteration (using new values for the 
primary controls), the transfer-function matrix was identified by successive, prescribed perturbations of 
the higher harmonic controls. During the controller iterations the matrix T was then recursively 
updated using an exponential filter, which weighted more heavily the most recent events (ref. 18). 

4. Results 

4.1 Hub Loads Correlations 

4.1.1 Rotor Performance and Trim Conditions 

All predictions of the rotor trim conditions and the hub loads shown in this section were obtained 
by trimming the CAMRAD model of the SA349/2 using the aircraft forces and moments (ref. 9). The 
correlation results are presented for four level-flight conditions ranging in advance ratio from 0.14 to 
0.37. Right parameters for the four conditions are listed in table 2. Figure 3 shows the difference 
between predicted and measured rotor trim angles, plotted as a function of advance ratio. Figure 3a 
indicates that the 1/rev longitudinal flapping angle, ^\c, is predicted to within 1/2 deg.; however, the 
shaft angle, and as a result the tip-path-plane angle (Og-Pic)/ are both in error by approximately 2 deg. 
CAMRAD’s prediction is nose up with respect to the measured data. Figure 3b shows the difference 
between predicted and measured rotor blade pitch trim parameters. In general, the collective and 
longitudinal cyclic pitch angles are overpredicted by an average of 1 deg. Lateral cyclic pitch is 
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overpredicted by 1/2 deg. The rotor blade trim data were obtained from the measured stick positions. 
These measurements may not be very accurate, however, because of flexibility in the linkages between 
the rotor blades and the control sticks. 

4.1 .2 Effect of the Rotor Wake Model on Hub Loads Correlations 

Correlations of the hub forces and moments in the shaft axis system were performed for the four 
flight conditions listed in table 2. The results, plotted as a function of advance ratio, are shown for both 
the rotating and nonrotating reference frames. The measured data are compared with predictions made 
using both a uniform inflow and a nonuniform inflow (prescribed wake) wake model. 

The hub force correlation, shown in figs. 4-6, was performed for the x-component of the total 
rotating in-plame force, f^x, the total resolved forces in the nonrotating frame, F\ and Fy, and the total 
nonrotating frame out-of-plane force, F^. figure 4 shows fjx 1/rev as a function of advance ratio, and 
the corresponding mean Fx and Fy components. An estimate of the SA349/2 hub drag and the X- 
component of the rotor weight were subtracted from Fx(mcas) and ftx(mea8) for each condition, to be 
consistent with the predictions. The estimates, which are listed in table 2, were derived from, 

Fx = Fx(meas) +mh g sinag - 0.5 fh p (Vq cosag)2, (13) 

where fh = 0.2 is the estimated D/q for the hub, mh is the mass of the rotor hub, Og is the shaft angle 
and Vq the helicopter velocity. Figure 4a shows that the in-plane rotating frame force decreased with 
advance ratio up to p^O.35. In the nonrotating frame, fig. 4 b shows that Fx also decreased with 
advance ratio, indicating that the lift vector became increasingly aligned with the aircraft eg. At low 
speed (p^0.14), where computed propulsive force was approximately 500 N, Fx contributed 
significantly to aircraft propulsion. At high speed (p;=0.37), however, the computed propulsive force 
increased to 3500 N, while Fx went to zero. Hence Fx was not a significant factor at high speed. Fy, on 
the other hand, increased with advance ratio, reflecting an increase in tail rotor thrust required to 
balance the rotor torque. The correlation of 1/rev ftx (fig- 4a) revealed a constant 200 N offset between 
the theory and the data. In the nonrotating frame, this offset was manifested in mean Fx at low speed 
and in mean Fy at high speed. The shift in offset from one fixed-frame force component to the other 
indicates a shift in the phase of f^x* Using a nonuniform inflow wake model offered no clear 
improvement to the correlation of steady drag and side force, because the trim state of the helicopter 
was the same for both wake models for the full advance ratio sweep. 

Figure 5 shows the f^x 2/ rev and 4 /rev rotating frame forces, and the corresponding 3 /rev 
non rota ting frame components of Fx and Fy. The fjx 2 /rev was the most significant contribution to the 
3 /rev fixed-frame loads, increasing with advance ratio. In the nonrotating frame, the magnitude of Fx 
3 /rev (predicted) was approximately 30% of the propulsive force magnitude. Both Fx 3/rev and Fy 
3 /rev exhibited the same trend of increasing with advance ratio. The correlation of oscillatory forces 
improved significantly with the use of a nonuniform inflow wake model. However, even using this 
wake model, the accuracy of the oscillatory force predictions generally deteriorated with increasing 
velocity. 

The prediction of 3/rev Fz is shown in fig. 6. The mean rotor thrust is not shown, since it was not 
accurately measured. The 3/rev Fz 3/rev correlation exhibited a significant improvement with the use 
of a nonuniform inflow wake model, although the measured force increased more rapidly with advance 
ratio than the predicted force. Comparing fig. 6b with fig. 5b reveals that the amplitude of 3/rev Fz 
was significantly greater than the amplitudes of the 3/rev in-plane forces, Fx and Fy. 

The hub moment correlation, shown in figs. 7-9, was performed for the total rotating in-plane 
moment, m^x/ the corresponding resolved components in the nonrotating frame, Mx and My, and the total 
nonrotating frame out-of-plane moment, Mz- Figure 7 shows the total mtx I / rev as a function of advance 
ratio, and the corresponding mean Mx and My components. As with ftx 1/rev, mtx 1/rev decreases with 
advance ratio, again suggesting that the lift vector and aircraft eg were aligning as the velocity 
increased. The trends of the fixed-frame mean moments were similar to the trends of the mean forces; 
the pitching moment (My) decreased, and the roll moment (Mx) increased with advance ratio. In both 
the rotating and nonrotating frame, the trend of the data was well predicted by CAMRAD. The 
correlation of 1/rev m^x (fig* 7a) improved with increasing advance ratio. In the nonrotating frame, the 
prediction of Mx was very good. Most of the rotating frame error was manifested in My. Using a 
nonuniform inflow wake model offered no clear improvement to the correlation in fig. 7, again reflecting 
that trim is independent of the wake model. 
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Figure 8a shows the mtx 2/rev and 4/rev rotating frame moments; the corresponding 3/rev 
nonrotating frame components of Mx and My are shown in fig. 8b. As with the hub forces, the mtx 2/rev 
is the most significant contribution to the 3 /rev moments, and all components increase with advance 
ratio. In the fixed-frame, the magnitudes of the measured 3/rev Mx and My are nearly identical, 
although the predicted Mx is higher than My. The correlation improved significantly with the 
introduction of a nonuniform inflow wake model, which reduced the error by a factor of two or more for 
3/rev My. 

The prediction of Mz is shown in fig. 9. The rotor torque, Mz (fig. 9a), was well predicted for all 
advance ratios, indicating a close match with rotor power as well. The prediction was within 7% using 
a uniform inflow wake model, however, using a nonuniform inflow wake model increased the 
correlation to within 1%. As expected, the mean torque increased with advance ratio. Figure 9b shows 
that the Mz 3/rev correlation exhibited significant improvement with the use of a nonuniform inflow 
wake model, particularly at the higher advance ratios. 

4.2 Effect on the Hub Loads of Changing the Blade Mass Distribution and Chordwise eg Location 

An examination of the influence of changes in the blade eg location and blade mass distribution on 
the hub loads was performed, in an effort to better understand the relationship between specific blade 
parameters and the subsequent forces and moments at the rotor hub. The highest speed condition, V3106 
(ji=0.37, Cx/ 0=0.066), was chosen for this investigation. Rgures 4-9 showed that the hub loads were 
generally largest for this condition. All predictions were made using a nonuniform inflow, prescribed 
wake model. 

Two changes in the blade mass distribution were studied, as shown in fig. 10. First, the blade mass 
was decreased 30% from 0-64% radial station, and increased 30% from 64-100% radial station, yielding 
a 12% net decrease in the total blade mass (MASS OUTBOARD). Second, the blade mass was increased 
30% from 0-64% radial station, and decreased 30% from 64-100% radial station, yielding a 12% net 
increase in the total blade mass (MASS INBOARD). 

Two changes in the eg location were also studied. The chordwise blade eg location was shifted 
forward 0.5% (CG FORWARD) and aft 2.0% (CG AFT) from the nominal 24.7% chord position. The 
forward position corresponded to the furthest forward position at which the analysis would converge. 
The aft position corresponded to the aft blade eg location flight tested by the SA349/2 helicopter (ref. 
14) in the second flight test. 

4.2.1 Overall Effect of Blade Structural Changes in the Hub Loads 

The overall effects of the changes in blade mass distribution and eg location on the amplitudes of 
the 3/rev forces and moments are shown in fig. 11. Moving the mass outboard decreased all the 3/rev 
loads. Moving the mass inboard increased the in-plane and decreased the out-of-plane forces and 
moments. Shifting the blade eg location from forward to aft increased the in-plane, and decreased the 
out-of-plane loads. Interestingly, the effect of the blade parameter changes was similar in trend for 
both the forces and moments, suggesting that the physical mechanisms for these changes were the same 
for both types of loads. Hence, the subsequent investigation of the effect on the hub loads of altering 
the blade properties focused on the hub forces only. 

4.2.2 Procedure for Investigating the Effect of Blade Parametric Changes in the Hub Loads 

Determining how and why alterations in the blade mass distribution and chordwise eg location 

affected the hub forces was accomplished by, 

1. Separating ftx and ftz into aerodynamic and inertial terms 

2. Examining the inertial terms for an individual blade to 

a. determine which bending modal components were affected by the changes in blade 

properties 

b. determine which inertial terms were responsible for the changes in 

modal components 

c. identify the sources of change; either through the modal characteristics 

(vk, *nk/ Iqk) o** through the aerodynamic forcing function 

(l/ac)(Mqk(aero)) 

3. Examining the aerodynamic terms for an individual blade to 

a. determine how single-blade hub force components (fca# fra/ fna) were affected 

b. determine the sources of change in f^/ through bending modal components 

c. determine the sources of change in fca and fna/ through blade lift 
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Following the above process, the next three sections investigate the effect of changes in the blade mass 
distribution and chord wise eg location on the hub forces. 

4.2.3 Aerodynamic and Inertial Components of fw and f^ 

To better understand how the blade alterations caused changes in the hub loads, the total rotating 
frame forces were first separated into inertial and aerodynamic components. The 3/rev in-plane 
nonrotating frame forces were derived from the 2/rev and 4/rev rotating frame forces. However, fig. 5 
showed that the predicted 4/rev contribution was much smaller than the predicted 2/rev contribution 
for the flight condition under investigation. Therefore, examining the effect of the blade parameter 
changes on the 2/rev rotating frame force was sufficient for determining the overall effect on the 3/rev 
nonrotating frame force. In addition, since the magnitudes of both ffx and fjy are equal, examining only 
ftx in the rotating frame was sufficient. The inertial and aerodynamic components of the 3 /rev out-of- 
plane force, ffz, were examined as well. 

Figure 12 presents vector plots of the total rotating frame 2/rev ftx hub force, separated into 
aerodynamic and inertial components. The effect of changing the blade mass distribution on the in- 
plane 2/rev force is shown in fig. 12a. The total vectors (T) are plotted for the baseline, the MASS 
OUTBOARD, and the MASS INBOARD cases, and each total vector has corresponding inertial (I) and 
aerodynamic (A) components. Figure 12a shows that the magnitude of the aerodynanuc vector is 
relatively insensitive to changes in the blade mass distribution. However, the phase of the 
aerodynamic component, and both the magnitude and phase of the inertial vector change (particularly 
for the mass inboard case). The increase of 2/rev fjx as it\e blade mass distribution shifts from outboard 
to inboard, and total blade mass increases 24%, is due primarily to changes in the inertied forces. The 
effect of shifting the blade chordwise eg location on the in-plane 2/rev force is shown in fig. 12b. 
Shifting the eg location from forward to aft results in a proportional increase of both the inertial and 
aerodynamic vectors, with no associated phase shift. 

Hgure 13 shows the breakdown of the out-of-plane 3/rev force into inertial and aerodynamic 
components. The effect of changes in the blade mass distribution is plotted in fig. 13a, where evidently 
most of the change in the total vector resulted from changes in the inertial component. Interestingly, 
the magnitude of the inertial vector decreased for both the MASS INBOARD and MASS OUTBOARD 
cases, causing the total vectors to decrease compared to the baseline. The effect of shifting the 
chord wise eg location is more straightforward, as shown in fig. 13b. The out-of-plane inertial and 
aerodynamic components both increased as the blade eg shifted from forward to aft. However, because 
the phase changed, the amplitude of the total vector decreased. 

Although the breakdown of hub forces into aerodynamic and inertial components reveals that 
altering the blade mass distribution affects primarily the inertial loads, while shifting the eg location 
influences both the aerodynamic and inertial components, the question of how these changes influence 
the aerodynamic and inertial forces has not been answered. Further examination of the aerodynamic 
and inertial components is required. A schematic of the rotating frame hub forces , representing the 
relationships between the single-blade forces (fc, fr/ and fn) and the total rotating frame forces (ftx and 
ftz)/ is shown in fig. 14. The schematic begins with the total rotating frame forces, ftx and ftz, splits 
them into their respective inertial and aerodynamic components (as already discussed), then shows the 
corresponding aerodynamic and inertial terms of the single-blade forces (fei, fca/ frl/ ^ra/ ^nl/ fna) as given 
in eqs. 4-6. The solid lines in fig. 14 lead to terms which are significantly affected by the blade 
structural changes under investigation. The other terms are either negligible contributions to the total, 
or unaffected by the changes in blade properties. The inertial terms can be decomposed into bending 
modal components. The aerodynamic terms are functions of a, which depends on torsion modal 
components. The aerodynamic term fra/ however, is a function of both a and the bending modes. 

4.2.4 Decomposing the Inertial Components 

The inertial terms of ftx and ftz were examined first. As shown by the expressions for the radial 
and chordwise inertial forces, frf and {q\, in fig. 14, the total inertial in-plane force, ftxl/ can be written 
as a sum of the modal components of fn and fei- Likewise the total inertial out-of-plane force, ftzi/ can 
be represented by the modal components of the single-blade inertial force, f|^. The significant modal 
components contributing to 2/rev f^xi and 3/rev ftzi are shown in fig. 15. The total inertial force for 
each case is labeled *T to be consistent with fig. 12. Hgure 15a indicates first of all that only the first 
three bending modes, rigid flap, rigid lag and first elastic bending, were significant to the total 2/rev 
in-plane inertial force. Changes in the 2/rev in-plane inertial force, ftxl/ caused by altering the blade 
mass distribution were manifested primarily through changes in the lag component of the first elastic 
bending modal contribution (k:=3). Recall that since coupled bending mc^es are used in the CAMRAD 
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analysis, each elastic bending mode has both a flap and a lag deflection component Changes caused by 
shifting the eg location were manifested through changes in the rigid lag bending contribution (k=l). 
The out-of-plane inertial force, 3/rev was virtually independent of the rigid lag bending mode 
(k=l). Figure 15b indicates that only the rigid flap and first elastic bending modal components 
contributed significantly to the 3/rev out-of-plane inertial force (I). The changes in 3/rev ft2i arising 
from changes in the blade mass distribution were manifested primarily through the flap component of 
the first elastic bending modal contribution (k=3), demonstrated in fig. 15b. Shifting the eg location, 
however, affected both the rigid flap (k=2) and first elastic bending (k=3) mode contributions to 3/rev 
^tzi* Table 3, which lists the blade bending modal frequencies for the five cases, reveals that changing 
the blade mass distribution affected the bending modal frequencies, while shifting the eg location did 
not. The torsion modal frequencies were unaffected for all cases. 

The next step was to answer the following two questions: 1) why does changing the blade mass 
distribution affect the first elastic coupled bending mode contribution to the 2/rev in-plane inertial 
force and the 3/rev out-of-plane inertial force, and 2) why does shifting the blade eg location affect the 
rigid lag contribution to the 2/rev in-plane inertial force and the rigid flap and first elastic bending 
mode contributions to 3/rev out-of-plane inertial force? Figure 16 shows a flowchart of all terms 
contributing to ftxi and ft^j. The inertial loads are a product of the modal response, and the blade 
inertial constants, Sqijc and Sqfk. The significant modal components discussed in fig. 15 are further 
decomposed in the flowchart of fig. 16. As fig. 16 indicates, the inertial constants, SqUe and Sqfk, 
depend on the mode shapes, iiik and iifk, while the blade modal degrees of freedom, qk# are functions of 
the blade modal frequencies, vk, the blade modal inertias, Iqk/ and the aerodynamic forcing function, 

(1 / ac)Mqk(aero)* Also, since the aerodynamic forcing function is itself a function of the blade bending 
mode shapes, qik, and Tjfk, and the single-blade aerodynamic forces, fca and fna/ then the blade modal 
degrees of freedom are a function of these terms as well. Thus the fundamental elements of the inertial 
forces are either modal inertial components: Vk, T|ik, Tjfk/ and Iqk, or aerodynamic components: fca and fna- 
(The aerodynamic components are decomposed in section 5.2.5.) The changes noted in the modal 
components of the inertial forces shown in figs. 15a and b were then traced to the modal inertial or 
aerodynamic elements. 

Beginning with the effect of changes in blade mass distribution on the in-plane inertial force, f^xi/ 
we examined the terms of the 2/rev first elastic bending modal component (labeled in fig. 16). The 
blade lag inertia, Sqik/ a function of the blade mass and the lag component of the bending mode shape, 

T|lk, were affected by changes in the blade mass distribution. Rgure 17a shows that the curvature of the 
lag component of the first elastic bending mode shape, ri[3, increased slightly as mass shifted outboard 
and decreased slightly when mass moved inboard. The 2/rev blade bending degrees of freedom, qk# 
plotted in fig. 17b for the first elastic bending mode, show the effect of each blade structural change on 
the first elastic bending degree of freedom, qa. Figure 17b reveals that qa was particularly sensitive to 
the blade mass distribution. The first elastic bending component of the aerodynamic forcing function, 
(l/ac)Mq3(aero), ^^Iso sensitive to blade mass distribution, as shown in fig. 17c. Since fig. 12a 
showed that the aerodynamic force was not significantly affected by altering the blade mass 
distribution, then the effect on the first elastic bending component of the aerodynamic forcing function, 
(l/ac)Mq3(aero)/ niust have arisen through the lag component of the first elastic bending mode shape, 

t] 13. Examining the other terms influencing the first elastic bending degree of freedom, the modal 
frequency, V3, and modal inertia, Iqs, we see in tables 3 and 4 that V3 and Iq3 also changed as the blade 
mass distribution shifted. Thus changes in the 2/rev in-plane inertial force, fjx/ resulting from changing 
the blade mass distribution, were caused by the first elastic bending modal elements: V3, Iq3 and 1113. 

The effect on the out-of-plane loads of changing the blade mass distribution was manifested in the 
same way, except through Sqf3^ t|f3 and the 3/ rev component of qs, as indicated in fig. 16. Rgure 18a 
shows that the flap component of the first elastic bending mode shape, , was sensitive to changes in 
the blade mass distribution. This component was also of much larger amplitude than the lag component 
shown in fig. 17a; therefore, this bending mode is more of a flap mode than a lag mode. Figure 18b 
reveals that changes in mass distribution affect the first elastic bending degree of freedom, q3, the most. 
Figure 18c shows that the first elastic bending component of the 3/rev aerodynamic forcing function, 
(l/ac)Mq3(aero)/ was also affected by the changes in the blade mass distribution. Since fig. 13a showed 
the insensitivity of the aerodynamic forces to blade mass changes, then the changes in 3/rev 
(l/ac)Mq3(aero) have arisen through the flap component of the first elastic bending mode shape. 
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T|f3. As with the 2/rev in-plane inertial force^ the 3/rev out-of-plane inertial force was affected by 
changes in the blade mass distribution through the first elastic bending modal elements: V3, Iq3 and T]f3 

A similar method was used to track down the significant effects on the in-plane 2/rev inertial 
force, ftxi due to shifting the eg location. Figure 15a has already shown that most changes in the in- 
plane 2/ rev inertial force arose from changes in the rigid lag modal component. A review of tables 3 
and 4 shows that the rigid lag modal frequency, V|, and modal inertia, Iq|, were not sensitive to blade eg 
location. The rigid lag mode shape. Tin, was also essentially unchanged. Hence the effect on the rigid 
lag component arose through the lag component of the 2/ rev aerodynamic forcing function, 

(1 /ac)Mqi(aero)- Figures 19 and 20 show that the magnitudes of both the rigid lag degree of freedom, qi, 
and the rigid lag component of the forcing function, (l/ac)Mqi(aeroV increased as the blade eg moved 
from forward to aft of the baseline position. This corresponds to the increase in the rigid lag component 
of 2/rev f^xi seen in fig. 15a. The effect of shifting eg location on the 2/rev in-plane inertial force was 
manifested through the aerodynamic terms of the forcing function, 2/rev f^a and 2/rev fna* 

The effect on the out-of-plane 3/ rev inertial force, caused by shifting the eg location was also 

manifested through the aerodynamic forcing functions, (l/ac)Mqk(aero)* Figure 15b showed, however, 
that the rigid flap (k=2) and first elastic coupled bending (k=3) modal components of were 
primarily affected (unlike ftxi/ where only the rigid lag mode was altered). The effects on the rigid 
flap component of both the rigid flap degree of freedom, q2, and the rigid flap component of the 
aerodynamic forcing function, (1 /ac)Mq2(aero)/ sre plotted in fig. 20. The amplitudes of the rigid flap 
components of q2 and (1 /ac)Mq2(aero)/ increased as the eg shifted from forward to aft. Revisiting figs. 
18b and c reveals that shifting the eg location had the opposite effect on the first elastic bending degree 
of freedom, q$, and the first elastic bending component of the aerodynamic forcing function, 
(l/ac)Mq2(aero)- These effects correspond to the trends seen in fig. 15b of the effect shifting the 
chord wise eg location had on the the rigid flap and first elastic bending modal components. The effect 
of shifting the eg location on 3/ rev f^zj was manifested through the aerodynamic forces of of the forcing 
function, 3/ rev fca and 3/rev f^a- The sources of change in these aerodynamic components are discussed in 
the next section. 

4.2.5 Decomposing the Aerodynamic Components 

The aerodynamic terms were examined in a manner similar to the inertial terms. The 2 /rev in- 
plane aerodynamic force, fjxa/ shown in fig. 12, was separated into its chordwise and radial single- 
blade aerodynamic force comp>onents, fca and fra- In addition, f|3, the radial aerodynamic force, was 
decomposed into its modal components. A vector plot of ftxa/ showing the contributions of fca/ the 
chordwise aerodynamic force, and the modal components of frg, is shown for the five cases in fig. 21. 
Figure 21 shows that most of the change in the 2/rev in-plane aerodynamic force, ftxa/ due to changes in 
the mass distribution, arise from changes in the radial force, fra- The changes in fra/ in turn, arise from 
changes in the first elastic bending modal component (k=3). Changes from shifting the eg location arise 
primarily from changes in the chordwise force, fca- The out-of-plane aerodynamic force was not 
similarly decomposed, since fig. 14 shows that f^za is not a function of the bending modes. 

Figure 22 presents a flowchart of the aerodynamic terms of fig. 14, and their components. Only the 
terms which were significantly affected by the changes in mass distribution and eg location are shown, 
and these have been decomposed into their basic m^al elements. As fig. 21 showed, changes in blade 
mass distribution affected the radial aerodynamic force, fra/ primarily through the first elastic 
bending modal component (k=3). Thus, altering the blade mass distribution affected the 2/rev in-plane 
aerodynamic force in the same manner as it affected the 2/rev in-plane inertial force, namely through 
the first elastic bending modal elements. 

The proportional change in the chordwise aerodynamic force, fca/ caused by shifting the chordwise 
eg location can be traced, according to fig. 22, to changes in the 2/rev component of the lift Likewise, 
the changes in the single-blade normal force, fna can be traced to changes in the 3/rev component of the 
lift The lift is a function of a, which is affected primarily by changes in the first elastic torsion modal 
component of the elastic twist, 00. Figure 23 supports this claim. Shown are vector plots of the 2/rev 
and 3/ rev components of the single-blade aerodynamic force normal to the hub plane, fna/ and the force 
normal to the hub plane due only to the first elastic torsion mode, fn(pl)/ where fn(pl) is given by, 

Wpl) = ^*5p (i^)^RLJ2cqpi Tcos 2<|> ^ dr/R (14) 
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Figures 23a and b show that the difference vector fna - fn(pl) remains unchanged from case to case, 
demonstrating that the changes seen in lift, and hence in fna/ arise from changes in the first elastic 
torsion mode. Furthermore, figs. 23a and b show that the single-blade aerodynamic force normal to the 
hub plane, fna/ is primarily sensitive to the aft shift in blade eg location. 

4.3 Effect of Higher Harmonic Control on the Hub Shears 

The HHC algorithm described in section 3.3 was used to reduce the predicted hub shears for the 
highest-speed flight condition (p;=0.37, Cj/c^O.066). The controller, or optimization algorithm, found 
a local minimum of the quadratic performance function, at which J was reduced by 95% from the 
uncontrolled value. The square root of the performance function,“jX is roughly the resultant 3/ rev hub 
force. At this high speed, the HHC required a large amplitude of blade pitch control. The following 
control vector was us^, 

^con = l(0O)3c (®0)3s (®c)3c Webs (®s)3c Wsbsl/ ('i 1 ) 

0con = [ -1.59, -0.22, -1.27, 0.26, -0.13, -2.12]. 

The pilot controls were given by 

0p = 13.13 - 2.45cos'F+ 7.56sin'F 

and the rotating frame HHC inputs corresponding to the control vector were 

0IIHC = -1.69cos2'F+ 0.19sin2M^- 1.59cos3'F- 0.22sin3'F+ 0.424T+0.19sin4T 

Time history plots of the rotating frame pilot control inputs, and the superimposed HHC inputs, are 
shown in fig. 24. The maximum HHC control input was a positive 2.8 deg. occurring at 4^100 deg. 

A bar chart showing the effect of HHC on the amplitudes of the force and moment coefficients is 
presented in fig. 25. Clearly, HHC managed to reduce the hub forces significantly. The hub moments, 
however, increased with the introduction of HHC. Thus, a performance function minimizing the 
amplitudes of all six loading components, rather than the amplitude and phase of just the hub forces, 
may be more desirable for reducing helicopter vibration. In addition, this analysis neglected shaft 
motion, which, if included, could change the hub shears produced by the rotor. The higher harmonic 
control also increased the power required by 10%. This increase was primarily due to an increase in 
profile power. Higher harmonic control increased the net angles of attack primarily on the retreating 
side of the rotor disk, thus increasing the profile power. This suggests also including the power in the 
quadratic performance function. 

Vector plots of the aerodynamic and inertial contributions to the fixed-frame hub shears, shown in 
fig. 26, illustrate the mechanism of hub shear reduction. Figure 26 shows the hub forces without HHC, 
with HHC, and the difference between the two, which is the perturbation of the forces due to HHC. In 
addition, fig. 26 shows the individual inertial (I) and aerodynamic (A) components of each total (T) 
vector. The total perturbation vector is created by generating large aerodynamic and rigid flap inertial 
perturbations. Thus, the hub force reduction is manifested not simply by direct cancellation of a 
particular hub-force component through generation of an opposing aerodynamic force, but by phase 
cancellation of large inertial and aerodynamic force components of the same order of magnitude. The 
effectiveness of HHC, therefore, depends not only on variations in blade angles of attack to achieve 
appropriate changes in the aerodynamic forces, but also on blade inertial characteristics to achieve 
appropriate changes in the inertial forces. 

5. Conclusions/Discussion 

Three topics related to the hub loads of the SA349/2 aircraft were investigated in this paper. 

First, correlations were presented between measured and predicted hub forces and moments, to examine 
the effect of the wake model and to develop confidence in the analysis. The last two subjects were 
theoretical investigations using CAMRAD and CAMRAD/JA. The effect on the hub loads of altering 
the blade mass distribution and eg location was examined. Finally, the effect of using higher harmonic 
control to reduce hub shear forces was studied. 

Using the CAMRAD analysis to predict measured hub loads resulted in the following conclusions: 


12 



1 . The 1 / rev rotating frame and mean nonrotating frame forces and moments were well 
predicted by the CAMRAD analysis using either a uniform inflow or a nonuniform inflow wake model, 
reflecting the good prediction of rotor trimmed forces and moments. 

2. Predictions of the 2/rev and 4/rev rotating frame and corresponding 3/rev fixed-frame 
forces and moments are significantly improved with the use of a nonuniform inflow model in the 
analysis. However, despite the significant improvement in the predictions of the oscillatory loads 
when using a nonuniform inflow model, the predictions are in error by approximately 30%, particularly 
at high advance ratio. 

Changing the blade mass distribution by shifting weight alternately outboard and inboard, and 
shifting the eg location aft and forward of 24.7% chord, affected hub loads at high speed as follows: 

1. The in-plane nonrotating frame 3/ rev hub forces and moments both increased when shifting 
the mass from outboard to inboard (with a corresponding 24% increase in the blade mass), and when 
shifting the eg location from forward to aft. 

2. The overall out-of-plane nonrotating frame 3/ rev forces and moments decreased when 
shifting mass both outboard and inboard, and when shifting the eg from forward to aft. 

3. Changing the blade mass distribution affected only the inertial part of hub forces, 
through changes in the first elastic bending modal components: the blade bending mode shape, the 
blade modal frequency, and the modal inertia. 

4. Shifting the blade eg location affected both the inertial and aerodynamic components of 
the hub forces. Chamges in the aerodynamic component were traced to changes in the first elastic 
torsion modal component. Changes in the inertial component were traced to the the aerodynamic 
forcing function, which was also influenced by the first elastic torsion mode. 

Intrtxlucing an HHC control algorithm into the CAMRAD /J A analysis had the following effects on 
the predicted hub shear forces for high speed flight: 

1. Using HHC, CAMRAD/JA predicted a reduction in the total resultant 3/rev hub force of 
95%, primarily by introducing large but opposing aerodynamic and inertial components of the hub 
forces. The evidence that HHC affects inertial as well as aerodynamic hub forces suggests that blade 
design optimization and HHC techniques should be considered together for the most beneficial results. 

2. The rotor power increased by 10% when HHC was introduced. The change was due to an 
increase in profile power caused by increased angles of attack on the retreating side of the rotor disk. 

These results are particular to the SA349/2 aircraft. In addition, this study of the effects of blade 
structural changes and higher harmonic control on hub loads is purely theoretical (no measurements 
were made). However, the approach of this investigation, which shows how the calculated results of 
a large computer code can be examined in detail, is very general. This type of study serves to elucidate 
the physical mechanisms for changes in rotor loads caused by blade (or other) design changes, which 
can be useful both for modifying rotor designs and improving rotor analyses. 
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Table 1. SA349/2 AIRCRAFT CHARACTERISTICS 


Aircraft gross weight 

Main rotor 

number of blades 

radius 

blade chord 

solidity 

blade tip speed 

Tail rotor blade tip speed 

Engine type 


19,600 N 


3 

5.25 m 
0.35 m 
0.064 
212 rpm 

215 rpm 

Astazou XIV H rated at 550 kw 


Table 2. Flight Parameters for Selected Conditions. 


Flight Condition 

V3I01 

V3103 

V3105 

V3106 

Altitude (m) 

306.2 

276.7 

317.2 

342.2 

TAS (m/s) 

30.2 

56.4 

74.4 

79.0 

Load Factor 

1.0 

1.0 

1.0 

1.0 

P 

0.14 

0.26 

0.35 

0.37 

Ct/cj 

0.0662 

0.0662 

0.0628 

0.0655 

Helicopter weight (kg) 

1999.0 

1991.0 

1979.0 

1973.0 

Static Pressure (mbar) 

978.0 

981.3 

976.7 

973.8 

Static Temperature (°C) 

9.2 

10.0 

11.2 

11.4 

Estimated hub drag (N) 

36.0 

305.0 

511.0 

597.0 


Table 3. Condition V3106. Bending, Vk, and torsion, COk. mode frequencies 


Mode 

k 

Baseline 

Mass 

outboard 

Mass 

inboard 

CG aft 

CG forward 

Rigid lag 

1 

0.58 

0.55 

0.63 

0.58 

0.58 

Rigid flap 

2 

1.03 

1.02 

1.03 

1.03 

1.03 

1st elastic 
bending 

3 

2.81 

3.30 

2.43 

2.81 

2.81 

2nd elastic 
bending 

4 

4.76 

5.31 

4.27 

4.77 

4.77 

3rd elastic 
bending 

5 

5.24 

5.83 

4.85 

5.23 

5.24 

1st elastic 

torsion 

1 

4.13 

4.13 

4.13 

4.13 

4.13 

2nd clastic 
torsion 

2 

12.81 

12.81 

12.81 

12.81 

12.81 
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Table 4. Condition V3106. Generalized modal inertias, Iqjj 


Mode 

k 

Baseline 

outboard 

Mass 

inboard 

Mass 

CG aft 

CG forward 

Rigid lag 

1 

0.86 

1.04 

0.68 

0.86 

0.86 

Rigid flap 

2 

0.98 

1.15 

0.81 

0.98 

0.98 

1st elastic 

3 

0.66 

0.61 

0.52 

0.66 

0.66 


bending 



Figure 1. Hub force and moment diagram. 
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Figure 3. Rotor trim predictions, 

a. difference between predicted and measured rotor disk orientation, 

b. difference between predicted and measured rotor control angles. 
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Figure 4. In-plane hub force correlations^ 

a. total 1/rev rotating-frame force, ftx* 

b. total mean nonrotating-frame forces, Fx and Fy, 
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Figure 5. in-plane hub force correlations, 

a. total 2/rev and 4/rev rotating-frame force, ftx. 

b. total 3/ rev noiuotating-frame forces, Fx and Fy. 
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Figure 6. Out-of-plane hub force correlation, 3/rev F^. 




Figure 7. In-plane hub moment correlations, 

a. total 1/rev rotating-frame moment, mtx. 

b. total mean nonrotating-frame moments, Mx and My. 
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Figure 8. In-plane hub moment correlations, 

a. total 2/rev and 4/rev rotating-frame moment, mtx. 

b. total 3/rev nonrotating-frame moment, M\ and My. 
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Figure 9. Out-of-plane hub force correlation, 

a. total mean moment, Mz, 

b. total 3/rev moment, Mz. 



Figure 10. Blade mass distributions. 
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Figure 11. Influence of changing the blade mass distribution and chordwise eg location on 
nonrotating-frame hub loads (|i=0.37, Ct/<J=0.066), 

a) Hub forces: Fx, Fy, Fz, 

b) Hub moments: Mx, My, M7. 












Figure 12. Vector representation of the total (T) rotating-frame force, 2/rev fjx, and the inertial 
(I) and aerodynamic (A) components, 

a. Effect of changing blade mass distribution, 

b. Effect of shifting chordwise eg location. 




Figure 13. Vector representation of the total (T) rotating-frame force, 3/rev ftz, and the inertial 
(I) and aerodynamic (A) components, 

a. Effect of changing blade mass distribution, 

b. Effect of shifting chordwise eg location. 
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Figure 14. Schematic of the aerodynamic and inertial terms. 
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Figure 15. Bending mode contributions to the total rotating-frame inertial hub forces (I), 

a. 2/rev ftxi, 

b. 3/rev ftzi. 
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Figure 16. Schematic breakdown of the inertial forces, showing terms influenced by changing 
the blade mass distribution and chordwise eg location. 
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Figure 17. Effect of changing blade properties on the 2/ rev first elastic bending modal 
component: 

a. lag component of the first elastic bending mode shape, 1113 , 

b. first elastic bending degree of freedom, q 3 , 

c. first elastic bending component of the aerodynamic moment, Mq3(aero)* 
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Figure 18. Effect of changing blade properties on the 3/rev first elastic bending modal 
component; 

a. flap component of the first elastic bending mode shape, rifa 

b. first elastic bending degree of freedom, qa 

c. First elastic bending component of the aerod)mamic forcing function, Mq3(aero)- 
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Figure 19. Effect of changing blade properties on the 2/rev rigid lag modal component: 

a. rigid lag degree of freedom, qi, 

b. rigid lag component of the aerodynamic moment, Mqj(aero). 
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Figure 20. Effect of changing blade properties on the 3/rev rigid flap modal component: 

a. rigid flap degree of freedom, q 2 , 

b. rigid flap component of the aerodynamic moment, Mq2(aero)- 
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Figure 21. Single-blade aerodynamic force contributions to the total in-plane rotating-frame 
aerodynamic hub force (A), including the modal components of fra- 



DISTRIBUTION 

Figure 22. Schematic breakdown of the aerodynamic forces, showing terms influenced by 
changing the blade mass distribution and chordwise eg location. 
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Figure 23. Vector plot of the single-blade rotating-frame aerodynamic force, fza» the first 
elastic torsion component of fza, fz(pl)» and the difference vector between the two, 

a. 2/rev, 

b. 3/rev. 
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Figure 23. Concluded. 
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Figure 24. Time history of the commanded control input, and the superimposed higher 
harmonic control command. 
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Figure 25. Influence of higher harmonic control on the amplitudes of the 3/rev nonrotating- 
frame hub loads, 

a. hub forces, 

b. hub moments. 



PERTURB. DUE TO HHC 

WITHOUT HHC 

WITH HHC 


Figure 26. Influence of higher harmonic control on the amplitude and phase of the total (T) 
nonrotating-frame hub forces, and the inertial (I) and aerodynamic (A) components, 

a. 3/reyFx^ 

b. 3/revFY, 

c. 3/rev Fz. 
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